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ABSTRACT 

We use high angular resolution data, measured from visibility of sources at the longest base- 
line of 4500 m of the Australia Telescope Compact Array (ATCA), for the Australia Telescope 
20 GHz (AT20G) survey to obtain angular size information for > 94% of AT20G sources. We 
confirm the previous AT20G result that due to the high survey frequency of 20 GHz, the 
source population is strongly dominated by compact sources (79%). At 0.15 arcseconds an- 
gular resolution limit, we show a very strong correlation between the compact and extended 
sources with flat and steep-spectrum sources respectively. Thus, we provide a firm physical 
basis for the traditional spectral classification into flat and steep-spectrum sources to separate 
compact and extended sources. We find the cut-off of -0.46 to be optimum for spectral indices 
between 1 and 5 GHz and, hence, recommend the continued use of -0.5 for future studies. 

We study the effect of spectral curvature on redshift cut-off of compact AGNs using 
recently published redshift data. Using spectral indices at different frequencies, we correct for 
the redshift effect and produce restframe frequency spectra for compact sources for redshift up 
to ~ 5. We show that the flat spectra of most compact sources start to steepen at ^ 30 GHz. At 
higher frequencies, the spectra of both populations are steep so the use of spectral index does 
not separate the compact and extended source populations as well as in lower frequencies. 
We find that due to the spectral steepening, surveys of compact sources at higher frequencies 
(> 5 GHz) will have redshift cut-off due to spectral curvature but at lower frequencies, the 
surveys are not significantly affected by spectral curvature, thus, the evidence for a strong 
redshift cut-off in AGNs found in lower frequency surveys is a real cut-off and not a result of 
K-correction. 
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1 INTRODUCTION 

Strong evolution of the number density of Active Galactic Nu- 
clei (AGNs) with redshift has been discussed in many papers. 
Studies done in optical. X-ray and radio frequencies have iden- 
tified a "quasar epoch", where the number density rises sharply 
as a function of redshift from z = to redshift of about 2.5 
to 3. The space density has been shown to decl ine beyond the 
redshift of 3 in fl at-spec trum radio sources (eg. |Peacocklll985l : 
lDunlop& Peac ock" 1990; 'Shaver et al."l996"; ' Wall et alj|2005h , in 
X-ray se lected q uasars (Hasinger et al. 2005) and o ptically selected 
quasars dSchmidt et alj|l99ll: IVigotti et all |2()03|) . The details of 
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whether such an evolution is purely due to space density evolution 
or due to luminosity evolution is yet to be resolved. Similarly, the 
decline in the number density in steep-spectrum radio sources has 
also been di s cusse d in the literature (eg. Dunlop & Peacock ( 19 9^; 
Ijarvis et alj ( 1200 lb ) and has been studied in detail by iRigbv et al.l 
l l201lh . 



The decline in the number density of QSO population be- 
yond redshift ~ 3 has been termed "redshift cut-otT' and has 
been explored in some detail by many authors. Some explana- 
tions offered suggest lumi nosity dependent space density evolution 
jPunlo p & Peacock 1990) while others suggest luminosity evolu- 
tion (e.g. Boyle, Shanks, & Peterson 1988). The rate of decline in 
the number density has al so be en discussed i n deta il. For example, 
iDunlop & Peacockl ( Il990h and lVigotti et"ai] l l2003h suggest a slow 
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decline while lShaver et alj ( Il996h and lWall et al] i2005h find a rapid 
decline in the number density of flat-spectrum radio quasars at high 
reds hift, comparable to results from X-ray and optical work. How- 
ever, |J^vis, feRawlingslfcood) have argued that spectral curvature 
causes the observable volume at high redshift universe to be too 
small to find any luminous flat-spectrum sources, removing the ev- 
idence fo r the rapid decline in number density at high redshift sug- 
gested bv lShaveretai] l ll996l) . The spectral curvature of sources at 
different redshifts gives the impression of redshift cut-off by caus- 
ing the sources to fall outside of selection criteria either by making 
the flux density fall under the selection limit or the spectral index 
to fall outside the selection limit, in popula tions selected bas ed on 
their spectra (eg. flat-spectrum population). IWall et al.l l l2005h have 
found that the use of non-contemporaneous flux density data can 
cause the effect of steep spectral c urvature similar to th at seen by 
Ijarvis & Rawlingsl ( I2OO0I) (see also lde Zotti et alj jZOlOl) ). 

The abundance of quasars at redshifts 2.5 - 3 has implications 
for the epoch of formation of the first structures in the universe, 
which has important consequences for the formation and evolution 
of massive galaxies. This epoch also has importance in determin- 
ing whether AGNs can be significant sources of reionization of the 
universe. 

Radio frequency studies of the redshift cut-off have the ad- 
vantage of being free from dust obscuration. Such studies have 
been done in the past with populations of flat-spectrum sources, 
selected from low frequency (^ 5 GHz) source catalogues. Flat- 
spectrum sou rces have traditionally been associated with com- 
pact sources ( IWalllll977t)^ with the sp ectral transition placed vari- 
ably around -0.4 (eg. Shaver 1995; Wall et al. 2005*), -0.5 (eg. 
Peacock 1985; Saikia et al. 2001; Browne et al. 2003) or -0.6 (eg. 
De Breuck et alj|2000l) . [n this paper, we follow the convention 
of Sv oc v"' and use -0.5 as the flat-spectrum cut-off. The flat- 
spectrum arises from the cumulative effect of multiple, very com- 
pact, self absorbed components peaking at different frequencies. 
Thus, compact AGNs are considered to have a very flat-spectrum, 
often spanning decades of frequencies. 

The Au s tralia Telescope 20 GHz (AT20G) survey 
( iMurphv et all I2OI0I) is a complete sample of radio sources 
and is do minated by comp act QSOs, due to its high frequency 
selection jSadler et alj2006! ). These sources have spectra measured 
at one to three different frequencies (20 GHz, 8.6 GHz and 4.8 
GHz) and have counterparts at one to two low frequencies (~ 1 
GHz) . The three high frequencies were measured near simulta- 
neously as part of the survey to remove any effects of variability. 
ISadler et al.l ( I2OO6I) find the median variability index at 20 GHz 
to be 6.9 per cent with 5 per cent of sources varying by > 30 per 
cent in flux density on a 1-2 year time-scale. The AT20G sources 
also have an angular size measurement at sub arcsec resolution 
which cleanly separates the compact AGNs from their jets and 
extended radio-lobes. In this paper, we demonstrate the strong 
correlation between compact sources and flat-spectrum sources 
providing a firm physical basis for the use of spectral index for 
such a separation. We then look at the radio spectrum of these 
compact sources to investigate the effect of redshift cut-off. 

In section 2, we outline the source sample used for this in- 
vestigation. In section 3, we explain the use of visibility to ob- 
tain an angular size measurement. Section 4 discusses the mea- 
surement of spectra of compact sources. Section 5 discusses the 
properties of such compact sources. In section 6 we investigate 
the effect of spectral curvature and its effect on the redshift cut-off 
of QSOs. Section 7 is the summary of this work. ACDM cosmol- 



Table 1. Number of AT20G sources with follow-ups, cross-matches with 
low frequency catalogues and optical counterparts. 



Population Nuinber 



iotat Al zulj 




6km visibility measured 


5542 


Ar20G^ 


3795 


AT20G^ + 6km visibihty -1- NVSS/SUMSS 


3403 


AT20G + Optical ID 


3873 


AT20G + Optical ID -1- Redshift 


1460 


AT20G + Optical ID -1- Redshift -1- 




6km visibility + NVSS/SUMSS 


1377 


+ Optical ID -1- Redshift + 




6km visibihty + NVSS/SUMSS 


838 



AT20G^ is the subpopulation of AT20G with 8.6 and 4.8 GHz follow-up. 
Sources with 5 > —15° do not have 8.6 and 4.8 GHz follow-up. 



ogy with Hp = 71 kms"^Mpc"\ Q.^ = 0.27 and = 0.73 
dLarson et al]l20T ij) has been used for this paper. 



2 THE AT20G SAMPLE 

The AT20G is a blind extragalactic survey carried out at a high fre- 
quency of 20 GHz for all declinations of the southern sky using the 
Australia Telescope Compact Array (ATCA). It covers a total area 
of 20, 086 deg'^. It has a galactic latitude cut-off |&| = 1.5 degrees. 
The survey has a total of 5890 sources above the flux density limit 
of 40 mjy at 20 GHz. The AT20G is the largest blind survey done 
at such a high radio frequency. Most sources south of the decli- 
nation of -15 degrees have near-simultaneous follow-up observa- 
tions within a month at 4.8 and 8.6 GHz (Murphy et al. 2010). The 
sources not observed at 4.8 and 8.6 GHz were a result of schedul- 
ing logistics and can not introduce any spectral bi as. The cross- 
match of AT20G sources against the 1.4 GHz NVSS jCondon etal] 
'1998*) catalogue for sources north of 5 = —40 degrees and 0.843 
GHz SUMSS ( Mauch e t al. 2003) catalogue for sources south of 
5 = — 30 degrees provides between one to two spectral points near 
1 GHz resulting in flux densities at two to five frequencies for the 
sample. Although the lower frequency points are not simultaneous, 
there is relatively little variability at these frequencies. Despite the 
galactic cut-off, the AT20G catalogue still contains a small num- 
ber of galactic thermal sources which are extended sources of op- 
tically thin (flat-spectrum) free-free emission. These high latitude 
sources are easily seen in Ha surveys and have been identified and 
removed from the main population for this analysis of the extra- 
galactic sources. 

Interferometric visibilities on long spacings (henceforth: 6km 
visibility) for > 94% of the AT20G sources have been measured 
(see section 3). 3403 AT20G sources have four or more spectral 
points between 1 and 20 GHz as well as 6km visibility. Optical 
identifications for 3873 AT20G sources hav e been made and red - 
shifts of 1460 sources have been compiled jMahonv et alj|2oTlh . 
This work makes use of the subpopulation of 838 sources with mea- 
sured redshift, 6km visibility and flux densities measured at four to 
five frequencies. Table[T]provides a summary of the different sub- 
populations. 
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Figure 1. The distribution of all 6km visibilities in the AT20G. Unresolved 
sources have a visibihty of 1 marked by the sohd vertical line. 



3 THE MEASUREMENT OF SOURCE SIZE 

The AT20G follow-up was done with the hybrid configurations of 
the ATCA, with five antennae in a compact configuration and the 
sixth antenna (called the 6km antenna in this paper) at a baseline of 
~ 4500m from the location of the compact configuration. This al- 
lowed us to calculate a visibility amplitude for over 94% of AT20G 
sources (Table [T). The remaining sources do not have their 6km 
visibility due to inclement weather during observation or unavail- 
ability of the 6km antenna due to maintenance. Thus, the sources 
with the visibility data are free from any bias. 

We calculate the 6km visibility from the ratio of the scalar 
average of the five long baseline amplitudes to the scalar average 
of the 10 short baseline amplitudes for the AT20G sources at 20 
GHz (Chhetri et al. in prep). Although, the AT20G follow-up ob- 
servations were done in a hybrid mode, the ATCA when used with 
its 6km antenna is essentially a linear array. So, the 6km visibil- 
ity corresponds to a position angle on the sky which is dependent 
upon the hour angle of observation. Most AT20G sources have vis- 
ibilities measured at two or more widely separated hour angles. 
Fig. [T] shows the distribution of all AT20G 6km visibilities for all 
sources. The 6km visibility for a point sources is expected to be 1. 
The Fourier transform provides a direct relationship between visi- 
bility and angular size of a source. Such a source size relationship is 
model dependent (eg. Gaussian, double peaked etc.) but, none the 
less, provides an effective method of separating compact sources 
from extended sources based on their size. 

Since it is not possible for any source to have visibility > 1, 
we can use the distribution of visibilities > 1 as an empirical es- 
timate of the effect of noise on the visibility estimate. Visibilities 
found to be > 1 in Figs. [T] and [3] are due to the effects of random 
noise on 6km amplitudes of unresolved sources with a RMS scatter 
of 0.053 to a maximum of 1.14. Considering a symmetric noise re- 
lated distribution from the unresolved value of 1, we arrive at 0.86, 
where we separate the compact sources from the extended sources 
and only 0.5% of compact sources would be misclassified as ex- 
tended. Therefore, for the purpose of this paper, we define sources 
with 6km visibility > 0.86 to be compact sources. For a Gaussian 
source, a visibility of 0.86 at 4500 m corresponds to a size of 0.15 
arcsecond at 20 GHz. 

We find that 79% of the sources in our survey are compact 



Figure 2. Redshift distribution of 6km visibility for flat (filled circles, 
Oj * > —0.5) and steep-spectrum (open squares, af'^ —0.5) sources. 
Visibility models ai'e drawn for 0.5, 1, 2 and 4 kpc diameter Gaussian 
sources. 

AGNs which is a result of the high survey frequency. This is very 
different from low frequency surveys such as NVS S which has only 
25% of the flat-spectrum AGN source population jPe Breuck et al.l 
l2000h . 

In Fig. |2] we have plotted the visibility on the 6km antenna 
baselines against redshift. Models for sources with an assumed 
Gaussian distribution of brightness are shown for a range of lin- 
ear sizes (full width at half maximum). If the source brightness 
distribution is double peaked rather than Gaussian the curves have 
exactly the same form but the sizes are scaled down by a factor of 
1.5. In reality there is variation in source morphologies and a single 
visibility cannot distinguish between these structures so the lines in 
Fig.[2]correspond to a small range of possible sizes. However, this 
uncertainty resulting from unknown source morphology is a rela- 
tively small fraction of the known range in radio source sizes from 
cores in the central few hundred parsec to megaparsec scale lobes. 

One interesting variation in the interpretation of the visibil- 
ity occurs for sources with emission coming from multiple regions 
with very different linear sizes. For example sources which have 
a significant flux in an unresolved core (<0.15 arcsec in our case) 
and the rest of the flux in a resolved jet or lobe structure. In this 
case the visibility should be interpreted as the fraction of the flux 
in the unresolved component, and the curves in Fig. |2] indicate the 
scale size for which the core component is unresolved. 



4 PROPERTIES OF COMPACT SOURCES 

Over the years, the spectral index of a radio source has been the 
primary method for selecting particular class of sources such as 
compact steep spectrum (CSS) sources, gigahertz peaked spectrum 
(GPS) sources etc. (eg. IO'Dej[T993 : ISaikia et alj|200lh : for se- 
lecti ng candidates for gravitationally lensed flat-spectrum sou rces 
(eg. iMvers et"ai]|l995l;1winn et al.|[200ol: ISrowne et alj|2003l) ; for 
the classification of populations such as beamed and unbeamed 
(eg. lackson & Wall 1999); f or modelling of source populations 
(eg. Dunlop & Peacock 1990: larvis & Rawhngsl l2000|- IWall et al.l 
[2005; de Zotti_eLal. 2005; Massardi et al. 2010); and for physical 
models (eg. lTucci et al.l20Tl]) . Fig.[3]plots spectral index between 1 
and 4.8 GHz against their 6km visibility for 3403 sources. With the 
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Table 2. Gaussian fits to the spectral index distribution of compact and ex- 
tended sources in Fig.|4] 



-1 1 

Spectral Index (1 to 5 GHz) 

Figure 3. Spectral index distribution plotted against 6km visibility. The dot- 
ted lines mark the limits of fiat v steep-spectrum on the x-axis and compact 
V extended sources on the y-axis. Solid Une at y = 1 represents the expected 
values of visibility for unresolved sources. The flux density for 1 GHz is 
obtained either from SUMSS or NVSS catalogue while 5 GHz is obtained 
from AT20G follow-up at 4.8 GHz. 



separation between compact and extended sources at the 6km visi- 
bility of 0.86, there is a very clean separation between compact flat- 
spectrum and extended steep-spectrum sources at the a ~ -0.5 line. 
The remarkably clean separation of compact flat-spectrum sources 
and extended steep-spectrum sources gives a very firm physical ba- 
sis for the traditionally used spectral index as a separator of com- 
pact and extended sources. Only ill out of 2347 (5%) flat spectrum 
sources are extended while 406 out of 3160 (12%) compact sources 
are steep-spectrum sources. 

Fig- El shows the distribution of spectral index {a\'^) for both 
compact and extended source populations. The mean and standard 
deviation for Gaussian fits for the compact and extended source 
populations are given in Table |2] 

We can use the spectral index distributions for compact and 
extended sources plotted in Fig.|4]to find the optimum value of the 
spectral cut-off to classify the flat and steep-spectrum AGN pop- 
ulations. We get slightly different answers depending on what we 
want to optimise. The maximum value of the joint probability that 
both flat and steep-spectrum sources are correctly classified is at 
a = -0.46 and 80% of the sources are correctly classified (see Fig. 
|4ll. This is independent of the relative number of flat and steep- 
spectrum sources in the survey, and consequently is only weakly 
dependent on the survey frequency. This estimate is right between 
the two values (-0.4 and -0.5) mostly used historically, for con- 
sistency we recommend the use of a -0.5 cut-off in future studies. 
However, if we want to maximise the total number of correctly clas- 
sified compact sources this does depend on the relative amplitude 
of the two populations and is dependent on survey frequency vary- 
ing between a = -0.4 for a survey frequency of 5 GHz to a = -0.6 
for a survey frequency of 20 GHz. 



5 THE SPECTRA OF SOURCES WITH 6KM VISIBILITY 

The overall spectral properties of AT20G sources sh ow com- 
plex s pectral changes and flux density dependence (Mass ardi et al.l 
l201lh . In Fig.|5]we compare the spectra for compact and extended 
sources for three different frequency ranges. We note the different 



Survey Frequency 


178 MHz 


20 GHz 


Data 


3C* 


AT20G Extended AT20G Compact 


Mean 


-0.85 


-0.75 0.00 


Sigma 


0.20 


0.23 0.42 



*3C values for ga laxies for spectral index between 0.75 and 5 GHz from 
iKellermann et al.l fT96ft) . 
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Figure 4. Distribution of Compact and Extended sources over spectral in- 
dex. The spectral indices are calculated between SUMSS/NVSS 1 GHz and 
AT20G 4.8 GHz flux densities. Mean spectral indices and dispersions from 
Gaussian fits to extended and compact source population distributions are 
given in Tabled] The solid curve is the joint probability of correctness (see 
text). The histogram for extended sources are shifted by -0.05 to show both 
populations at the steep-spectrum end. 



characteristics between the three frequencies. Spectral index be- 
tween two frequencies with larger separation is less susceptible to 
errors due to flux measurement errors than spectral index between 
two frequencies close to each other For a 10% error in flux density 
in Fig.|5] the af^ has an error of 0. 16, the ag% has an error of 0. 1 1 
and af has an error of 0.06. The 5 and 8 GHz flux density mea- 
surements were done within a month of the 20 GHz measurements. 
Hence, they do not suf fer much from variability of compact sources 
dMassardi et al.ll201ll) . Although the 1 GHz flux densities were 
measured at different epochs (between minimum of one to max- 
imum of 15 years), variability is much less at lower frequencies. 
Hence, errors in flux densities or due to variability do not explain 
most of the differences seen in the three plots. 

In Fig. [5] we note the appearance of a steep-spectrum 'tail' 
for compact sources in the top a|"g plot. This indicates that while 
compact sources are mostly flat spectrum at lower frequencies, a 
fraction of compact sources have steeper spectra between 8 and 
20 GHz. It is also evident that the population of compact inverted 
spectrum sources, sources with spectral index > 1, at lower fre- 
quencies shifts to lower spectral index and combines with the main 
population of the compact sources at higher frequencies. These are 
compact core sources that are rising at lower frequencies and have 
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Figure 6. Spectra of compact and extended sources with redshift calcu- 
lated between 8.6 and 20 GHz. Filled circles represent compact sources and 
open triangles represent extended sources. The median and standard error 
for redshift bins, marked by dashed veilical lines, are shown in Table|3] 



Figure 5. 6km visibility is plotted against three spectral indices increasing 
in frequencies with a*-* at the bottom panel, o| | at the mid panel and 
Og'g at the top panel. The solid vertical line corresponds to spectral index 
= -0.5 and the solid horizontal lines correspond to 6km visibility = 1. 

either become flat or steep at higher frequencies (e.g. GPS sources) 
dSadleret aljioog) . 

On the other hand the spectra of extended sources which are 
well separated by spectral index at low frequencies, have a larger 
scatter at higher frequencies. A significant fraction of the extended 
sources move towards the flatter part of the plot at higher frequency. 
These extended sources are mostly composed of a flat-spectrum 
core with a steep-spectrum jet, lobe or a hotspot. With increasing 
frequency, the steep-spectrum part of the sources gets weaker while 
the compact flat-spectrum cores start to dominate the spectrum. 
These composite sources show an overall combined flat spectrum 
but are still extended. This domination of flat-spectrum core should 
be easily observable by imaging these composite sources at higher 
frequencies, when only the compact cores should be visible. With 
increasing frequency, there is a further steepening in the spectra of 
the extended structures due to energy loss. This is evident in the 
left part of Fig. |5] with the appearance of a population of very steep 
extended sources in the q|'^ and plots. These two effects com- 
bined produce the larger scatter seen in the a||g plot of Fig. [5] for 
extended sources. 

It should also be noted that at higher frequencies the spectra 
of compact sources starts to become indistinguishable from steep- 
spectrum sources (Fig. |5] al'5 plot). It is then clear that while the 
use of spectral index to separate compact and extended sources is 
valid using lower frequency spectra, it is unable to provide a clean 
separation of compact and extended sources at higher frequencies. 



6 THE REDSHIFT CUT-OFF AND THE EFFECT OF 
SPECTRAL CURVATURE 

IWall et alj j2005l) used 379 radio selected QSOs from the Parkes 
quarter-Jansky flat-spectrum sample selected at 2.7 GHz to study 
the evolution of flat-spectrum quasars at high redshift. They found 



that the redshift cut-off is real and is similar to the cut-off observed 
in optically selected QSOs. Similar res ults have been reported fo r 
cosmic star formation rates. However, Ijarvis & Rawlingsl ( I2000I) . 
using the Parkes half-Jansky flat-spectrum sample selected at 2.7 
GHz, argue that because of spectral curvature and the resulting K- 
correction there is not enough observable volume at high redshift 
to find the brightest flat-spectrum radio sources, which gives the 
impression of redshift cut-off. We can use the sub population of 
compact sources from the AT20G survey, with redshifts, to study 
the effect of spectral curvature and the K-correction on the redshift 
cut-off. 

A survey at any frequency is biased to the sources that are 
stronger at the survey frequency. Since this is very obvious, it 
is not usually considered a bias. At lower frequencies the steep- 
spectrum sources dominate the population at all luminosities and 
only contain the rare extremely luminous flat-spectrum sources. 
Since it is the flat spectrum population that has a high quasar iden- 
tification rate, this is the sample used to investigate evolution of 
flat-spectrum quasa rs (Du nlop & Peacock 1990: Shaver et alj 19961; 
IJarvis & Rawlingdl2000l : rvigotti et alj|2003l : IWaU et alj|2005l) . Our 
frequency selection provides a much larger fraction of these sources 
without introducing a different class of sources, thus, w e reduce 
the bias by having a high frequency selected sample (e.g. ICondonI 
I2OOM . 

As noted in Section 2, 838 out of 1460 of AT20G sources with 
redshift have 4.8 and 8.6 GHz flux densities, 1 GHz NVSS/SUMSS 
counterparts and 20 GHz flux densities as well as 6km visibility. 
We use this subsample of sources for the following analysis. Of 
these, 683 sources are compact by our definition while the remain- 
ing 155 sources are extended. The distribution of the redshifts for 
this AT20G subpopulation is listed in Table[3] 

Fig. [6] plots the spectral indices for compact and extended 
sources between 8.6 and 20 GHz against redshift. It is evident from 
the plot that extended sources do not show apparent changes in the 
spectral index over the redshift range. Compact sources, however, 
show a distinct steepening of their spectra with redshift. It is appar- 
ent that a spectral classification is no longer separating the compact 
and extended sources above redshifts ^ 1 . 

Since most extended radio sources have close to power-law 
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Figure 7. Rest frame spectral index of compact and extended sources. Filled 
circles represent compact sources and open triangles represent extended 
sources. All points have been brought close to the rest frame spectral in- 
dex between 8.6 and 20 GHz by using the spectral index between 8.6 and 
20 GHz for all sources below z = 1, using spectral indices between 4.8 and 
8.6 GHz for all sources between redshifts 1 and 2 and using spectral index 
between 1 and 4.8 GHz for all sources with z ^ 2. 



spectra over a large range of radio frequencies, we do not expect 
a significant change in spectrum with redshift due to the change in 
restframe frequency. But as discussed in section 5, compact sources 
have flat spectra at lower frequencies that gets steeper at higher 
frequencies. Hence, the spectral indices of compact sources are 
affected b y the redshift effec t causing the redshift cut-off as dis- 
cussed by Ijarvis & Rawlingsl llOOdl . To test this idea further we 
use the spectra at different observed frequencies for different red- 
shift ranges. 

The median spectral indices and errors for these redshift 
ranges are hsted in Table [3] Our sources have flux densities mea- 
sured between four to five frequencies which allows us to calcu- 
late, at least, three contiguous spectral indices {oi\'^ , at;!, al^e)- 
We take advantage of these spectral points by using the best rep- 
resentative spectra for different redshift bins to remove the effect 
of redshift and, thus, compare their spectral indices between corre- 
sponding rest frequencies. For this, we approximate that the spec- 
tral index of af best represents the redshift bin between 2 and 5 
(rest frequencies ~8 and ~14 GHz), a| ;| best represents the red- 
shift bin between 1 and 2 (rest frequencies ~10 and ~26 GHz) and 
Qg^g represents the redshift bin for redshifts up to 1. Plotted in Fig. 
[T] this essentially brings all spectral indices to Qf|"g for compari- 
son in their restframes. We use this simple procedure to preserve 
the measured values rather than fitting complex spectral shapes to 
a small number of frequencies. 

As is evident in Fig. |7] the redshift correction completely re- 
moves the steepening of the spectral indices at higher redshift seen 
in Fig. [6] Table |3] even indicates that there is now an overall flat- 
tening of spectra at higher redshift to 9 sigma confidence. These 
results show that the spectral curvature plays a dominant role at 
higher frequencies. To illustrate this further, we have constructed 
the K-corrected spectra for the entire sample of sources in Fig. [8] 
Of the sub population of 1377 sources (see Table[Tll, 1092 sources 
are compact by our definition and only 671 sources have spectra 
at four to five frequencies. We calculated the rest frequencies for 
these sources and normalized them with their average flux density. 



and plotted them against the rest frequency in Fig. [8] The mean red- 
shift for these 671 co mpact sources is ^ 1. The figure confirms the 
results seen in fig. 6 o f lWalletalJ ( l2005l) which shows that the spec- 
tra of flat-spectrum sources at their rest frames at lower frequencies 
is predominantly flat suggesting that spectral curvature does not 
play a significant role at lower frequencies. At lower frequencies 
up to 5 GHz, we observe that there is a mixture of steep, flat 
and inverted spectrum sources. It is also clear from the figure that 
the majority of compact sources remain flat-spectrum up to ~30 
GHz in their rest frame frequencies at all redshifts but the trend is 
then for almost all sources to get steeper. The median normalized 
flux density for the rest frequency bin between 35 - 45 GHz drops 
by ^^15% compared to that for the rest frequency bin between 15 
- 25 GHz. This explains the effect of spectral steepening at high 
frequency of compact, and hence, flat-spectrum sources at higher 
redshifts seen in Fig. |6] and Table [3] as well as the removal of the 
effect after redshift correction in Fig. [7] The plot also indicates that 
at lower frequencies, the spectral curvature does not play a signif- 
icant role. In fact, there is a population of steep-spectrum core-jet 
sources at the lower frequency end of Fig. [8] which produces the 
opposite effect. As discussed in Section 5, these compact sources 
that have steep-spectrum jets at low frequencies are increasingly 
dominated by their flat-spectrum cores at higher frequencies. The 
redshift effect of these sources, then, makes them appear flat spec- 
trum at observed lower frequencies at higher redshifts. For sources 
at high redshifts, this effect negates the effect of spectral steepening 
of the overall spectra, at lower frequencies. 

Thus, the redshift cut-off is affected by the spectral curvature 
of compact sources as suggested by Jarvis et al., (2001), but this 
effect only becomes important at higher frequencies (f >5 GHz) 
and can be avoided by measuring the spectra at lower frequencies. 
Fortunate ly, this is where majority of evolution stud i es have been 
done ( e g.lPeacoc k (1985); Dunlop & Peacock ( 1990):l shaver et al.l 
( Il999h ; |jarvi7& Rawlin gs ( 2000,) ; .Wan et al. (2005)). 

So far, we have only discussed the effect of the spectral curva- 
ture on the redshift cut-off caused by the use of the spectral index 
to filter out the compact flat-spectrum sources. However, it is also 
necessary to see how this effects the nu mber of sources above the 
flux limit of the survey, referred to bv iJarvis & Rawlingsl ( |2000|) 
as the observable volume. This cut-off is determined by the survey 
frequency regardless of the frequency used to measure the spec- 
trum. Low frequency surveys will select against the compact flat- 
spectrum sources at all redshifts. Higher frequency surveys will 
have much larger populations of compact sources but they will also 
have a much stronger redshift cut-off due to spectral curvature. 

To evaluate this survey sensitivit y effect, w e plot the average 
parameters used in the i Jarvis & Rawlingsl ( I2OOOI) models, normal- 
ized at 2.7 GHz, with our rest frame spectra in Fig.[8] The model is 
given by a polynomial of the form y = logioSi, = E^^Qii^;', where 
X = logio(!//GHz) with ai = 0.07 ± 0.08 and 02 = -0.29 ± 0.06. 

It is apparent that the model predicts a much stronger spectral 
curvature than that obtained for the AT20G compact sources. Even 
the closest fit, at the extreme of the parameters range, is steeper 
than that obtained for the AT20G compact sources. We note that the 
data available to Jarvis & Rawlings (2000) was strongly biased to 
unusual spectra and GPS sources by the spectral data available and 
is not representative of the tru e compact source populati on. On the 
basis of their model spectrum I Jarvis & Rawlingsl ( I200CI) predicted 
a sharp cut-off in the number of observable sources at h igh red- 
shift a nd consequently argued that the cut-off seen bv Shaver et al.l 
d 19961) was not evidence for cosmological evol ution. On average, 
our spectra are at least five times brighter than Ijarvis & Rawlingsl 
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Table 3. Median spectral indices and standard error for compact and extended sources before (Fig. |6) and after (Fig.|7) redshift corrections. There are only 
two extended sources in the z ^ 2 bin. The median and standard error are kept inside parentheses to note their limited significance. 



Redshift 
Bins 




Compact 








Extended 






Total 


Number 


Uncorrected Error 


Corrected 


Error 


Number 


Uncorrected 


Error 


Corrected 


Error 


z < 1 


363 


-0.17 0.02 


-0.17 


0.02 


138 


-0.79 


0.05 


-0.79 


0.05 


501 


1 ^ z <2 


211 


-0.28 0.03 


-0.07 


0.03 


15 


-0.80 


0.08 


-0.73 


0.07 


226 


z > 2 


109 


-0.49 0.04 


0.10 


0.05 


2 


(-0.57) 


(0.38) 


(-0.24) 


(0.02) 


111 


Total 


683 








155 










838 



model, at our higher frequencies which are relevant for the 
high red shift sources found in lo wer frequency sampl es such as 
those of Ijarvis^ Rawlingd ilOOd) or lWall et alJ (l2005h . S)0, even 
at high redshift the expected number of sources will be larger by 
a similar fraction at this flux density level. The bias introduced by 
preferential selection of fiat-spectrum sources due to high AT20G 
survey frequency is appropriate for the high re dshift sources where 
cut-off is seen. Hence, the cut-off reported bv lshaver et alj ( fl996l) 
must be mainly due to cosmological evolution and is not the ef- 
fect of K-correction on either the number of sources or the abil- 
ity to make spec tral classification. The more detailed analysis of 
IWalletalJj200^ . which incorporated the various selection effects 
and biases explicitly, also demonstrated the reality of the cosmo- 
logical evolution component. 

We compared the Parkes quarter-Janksy flat-spectrum sam- 
ple ( |jacksonetalj |2002h with the AT20G data in the overlapping 
area of the sky. The PQJFS sample has 804 southern hemisphere 
sources and actually includes sources with flux densities less than 
250 mjy. We have included all sources down to 220 mjy since 
that corresponds to 100 mJy at 20 GHz for the spectral index 
cut-off of -0.4. There are 716 PQJFSS sources over the limit of 
220 mJy at 2.7 GHz. The cros s match of thes e sources' positions 
against the AT20G catalogue jMurDhvetal] |2010) finds 640 of 
these sources and further 3 1 are fou nd in the AT20G scanning sur- 
vey catalogue llHancock et al.ll201lh which avoids the incomplete- 
ness introduced by the AT20G follow-up obs ervations. We used th e 
updated positions of the sources provided in' J ackson et alj j2002h . 
where available, since the original PQJFSS positions are not fully 
reliable. Thus, our cross-match found a total of 671 (94 %) PQJFSS 
sources in the AT20G survey. Of the rema ining 45 sour c es, 14 
were later found to be extended sources by Ijackson et alj ( l2002h 
and are not included in their analysis. For the remainder, our lit- 
erature search finds that 7 are detected at lower frequencies and 
have steep-spectrum, 9 are known variable sources based on mul- 
tiple epoch observations and 1 source that has not been detected in 
any other surveys. The remaining 14 sources (2%) are well within 
the completeness expected from long term variability noting that 
the AT20G survey was done > 20 years after the PQJFS survey. 
Only 3 (0.4%) of these sources have evidence for extreme spec- 
tral steepening at high frequency. Hence, we can conclude that we 
have not lost a significant fraction of sources as a result of spectral 
steepening which has put them below our survey limit. 

603 of the detected sources have 6km visibilities and 548 of 
them are compact while 55 (9%) are extended, consistent with the 
overlap in the spectral index distribution shown in Fig.|4] Out of the 
detected sources, 461 sources have redshift information. The pub- 
lished PQJFS sample has 459 sources in the southern hemisphere 
with redshifts. 




Rest Frequency (GHz) 

Figure 8. Rest frame radio spectra of 671 compact AT20G sources. Sources 
have flux density observations at a subset of 0.8, 1.4 and 4.8, 8.6 and 20 
GHz. All observed frequencies are converted to the rest frame and the flux 
densities are normalized by their average flux density over all frequencies 
present. The solid black line is the fit to the spectra of the most lumi- 
nous P arkes Half-Jansky Flat-Spectrum sample done bv ljarvis & Rawlingj 
i2OO0l) normalized at 2.7 GHz survey frequency. 



7 SUMMARY 

We summarise the results of this work. 

(i) 6km visibility, calculated as the ratio of five long baselines to 
the 10 short baselines of the ATCA in hybrid configurations, pro- 
vides information on the angular sizes of the extragalactic source 
population at 20 GHz. It is used to separate the source population 
into compact and extended sources at subarcsec (0.15 arcsecond) 
level. 

(ii) 6km visibility provides a firm physical basis for spectral 
classification into flat and steep-spectrum sources to select compact 
and extended sources. For spectral indices measured between 1 and 
5 GHz, the spectral index of -0.46 is the optimum cut-off between 
the two populations. However, the traditional cut-off of either -0.4 
or -0.5 are almost as good and for consistency, we recommend that 
-0.5 be used in future studies. 

(iii) 79% of the sources in our survey are compact AGNs which 
is a result of the high survey frequency. This is very different from 
low frequency surveys such as NVSS which has only 25% of the 
flat-spectrum AGN source population. 

(iv) A larger scatter in the spectral index of extended sources is 
observed at higher frequencies due to the domination of compact 
cores in core-jet sources as well as the steepening of extended jets 
due to energy loss. Hence, the use of spectral index at high fre- 
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quencies is not as effective to select compact and extended sources. 
Using the radio spectrum below a few gigaliertz avoids this prob- 
lem. 

(v) The spectral curvature is significant for compact sources, but 
it only becomes significant at higher frequencies. Below a few giga- 
hertz frequencies, the redshift effect on compact sources is not sig- 
nificant. We note that the general spectral steepening seen above 20 
GHz will occur at observed frequencies around 5 GHz, for sources 
with redshift > 3. Hence, surveys at frequencies > 5 GHz will have 
a redshift cut-off due to spectral curvature. 

(vi) At high frequencies, the spectra of both compact and ex- 
tended sources are steep and the two populations become indistin- 
guishable from each other due to the spectral steepening of compact 
sources by ~ 15% between 20 and 40 GHz. 

(vii) The effect of K-correction shown by Ijarvis & Rawling^ 
2000) is i mportant at high fr equencies as discussed bv lShaver et al 



19960 an d lWall et al.l ( l2005h for flat-spectr um sources and also for 
steep-spectrum sources teigbv et al.l l201 1() . Ho wever, this should 
not affect the results of Shaver et al.l( ll996h and lWall et all j2005l) 
as they used low frequency spectral index to select compact 
sources. Hence, the conclusion that the compact sources show a 
strong redshift cut-off is supported by our data. 

In future, with the information on the completeness of the optical 
identification and redshifts, it will be possible to explore the density 
evolution and redshift cut-off of compact AGN without using the 
spectral index to select compact sources. 
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